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Abstract. We study a multiobjective variational problem on time scales. 
For this problem, necessary and sufficient conditions for weak local Pareto 
optimality are given. We also prove a necessary optimality condition for the 
isoperimetric problem with multiple constraints on time scales. 

1. Introduction 

The calculus on time scales was initiated by Aulbach and Hilger (see e.g. [2]) in 
order to create a theory that can unify discrete and continuous analysis. Since then, 
much active research has been observed all over the world (see e.g. [TJ |31 IH [7] and 
references therein). In this paper we consider multiobjective variational problems 
on time scales (Section l3.2p . By developing a theory for multiobjective optimization 
problems on a time scale, one obtains more general results that can be applied to 
discrete, continuous or hybrid domains. To the best of the authors' knowledge, no 
study has been done in this field for time scales. The main results of the paper 
provide methods for identifying weak locally Pareto optimal solutions; versions for 
continuous domain one can find e.g. in [HllHlll]- We show that necessary optimality 
conditions for isoperimetric problems are also necessary for local Pareto optimality 
for a multiobjective variational problem on a time scale (Theorem 13. 8|) . and the 
sufficient condition for local Pareto optimality can be reduced to the sufficient 
optimal condition for a basic problem of the calculus of variations on a time scale 
fTheorem l3.7p . We also prove a necessary optimality condition for the isoperimetric 
problem with multiple constraints on time scales (Section l3.ip . 

2. Time scales calculus 

In this section we introduce basic definitions and results that will be needed for 
the rest of the paper. For a more general theory of calculus on time scales, we refer 
the reader to [S] . 

A nonempty closed subset of K is called a time scale and it is denoted by T. 
The forward jump operator cr : T — > T is defined by 

a{t) = inf {s eT ■.s>t}, for aU t e T, 

while the backward jump operator p : T ^ T is defined by 

p{t) = sup {s e T : s < t}, for aU t e T, 

with inf = supT (i.e. a{M) = M if T has a maximum A/) and sup0 = inf T (i.e. 
p{ni) = m if T has a minimum m). 
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A point t G T is called right-dense^ right- scattered^ left-dense and left- scattered 
if a{t) = t, a{t) > t, p{t) = t and pit) < t, respectively. 

Throughout the paper we let T = [a, 6] fl Tq with a < b and Tq a time scale 
containing a and b. 

Remark 2.1. The time scales T considered in this work have a maximum b and, by 
definition, a{b) = b. 

The graininess function : T — s- [0, oo) is defined by 
p{t) = cr{t) - t, for aU t £ T. 

Following F, we define T*-' = T\{p{b),b], T'^'" = (T'^)''. 

We say that a function / ; T ^ M is delta differ entiable at f G T'^ if there 
exists a number f^{t) such that for all £ > there is a neighborhood [/ of < (i.e. 
U ={t-5,t + 5)r\T ioY some 5 > Q) such that 

\f{a{t)) ~ f{s) - f''{t){a{t) - s)\ < e\a{t) - s\, for aU s e U. 

We call (t) the delta derivative of / at t and say that / is delta differentiable on 
T*^ provided f^{t) exists for all teT''. 

For delta differentiable functions / and g, the next formula holds: 



where we abbreviate here and throughout the text f o a hy . 

A function / : T — > M is called rd- continuous if it is continuous at right-dense 
points and if its left-sided limit exists at left-dense points. We denote the set of all 
rd-continuous functions by Crd and the set of all delta differentiable functions with 
rd-continuous derivative by C^^^j. 

It is known that rd-continuous functions possess an antiderivative, i.e. there 
exists a function F with — /, and in this case the delta integral is defined 



by f{t)At ^ F{c) - F{d) for all c, d e T. The deha integral has the following 



We now present the integration by parts formulas for the delta integral: 
Lemma 2.2. (^) Ifc,dGT and f,g GC],^, then 



We say that / : T — > R" is a rd-continuous (a delta differentiable) function if each 
component of /, : T — > M, is a rd-continuous (a delta differentiable) function. 
By abuse of notation, we continue to write Crd for the set of all rd-continuous 
vector valued functions and Cj^j for the set of all delta differentiable vector valued 
functions with rd-continuous derivative. 

The following Dubois-Reymond lemma for the calculus of variations on time 
scales will be useful for our purposes. 

Lemma 2.3. (Lemma of Dubois-Reymond ^) Let g E Crd, 9 ■ [a, b]'' M". Then, 



if9)^it)^f^{t)g^t) + f{t)g^{t) 
^f^it)9{t) + r{t)g^{t), 



property: 






git) ■ r]^{t)At = for all rj G C].^ with 77(a) = 77(6) 
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if and only if g{t) — c on [a, 5]*^ for some c € M". 

3. Main Results 

We begin by proving necessary optimality conditions for isoperimetric problems 
on time scales (i j3.ip . In t j3.2l we show that Pareto solutions of multiobjective vari- 
ational problems on time scales are minimizers of a certain family of isoperimetric 
problems on time scales. 

3.1. Isoperimetric problem on time scales. 
Definition 3.1. For / : [a, b] — > M" we define the norm 

niax ||r(0||+ niax 11/^(011, 

where || • || stands for any norm in R". 

Let C : C^^ ^ M be a functional defined on the function space C^^ endowed with 
the norm || ■ \\(^a^ and let A C C^^. 

Definition 3.2. A function / G A is called a weak local minimum of C provided 
there exists S>0 such that £[/] < £[/] for all f £ A with ||/ - f\\ci^ < 6. 

Now, let us consider a functional of the form 
(1) C[y]^ f L{t,y%t),y^{t))M, 

J a 

where a, 6 G T with a < L{t^ s, v) : has partial continuous 

derivatives with respect to the second and third variables for all t G [a, 6]*^, and 
L{t, •, •) and its partial derivatives are rd-continuous at t. The isoperimetric problem 
consists of finding a function y satisfying: 

(1) the boundary conditions 

(2) y{a)=a, y{b)^p, a,/3eR"; 
and 

(ii) constraints of the form 

(3) gM= G,{t,y^{t),y^{t))At = (,, i = l,...m, 



where £,i,i = 1, ... to, are specified real constrains, Gi{t, s, v) : [a, b]'^ xR" xR" — > R, 
i = 1,...TO, have partial continuous derivatives with respect to the second and 
third variables for all t G [a,b]'^, and Gi{t,-,-) and their partial derivatives are 
rd-continuous at t] that takes ([1]) to a minimum. 

Definition 3.3. Let £ be a functional defined on C^^. The first variation of C at 
y G C^^ in the direction rj G C^^, also called Gateaux derivative with respect to r] 
at y, is defined as 

SRy; m = hm^^o = -^^y + sri\\e=o 

e oe 

(provided it exists). If the limit exists for all rj G G^^, then C is said to be Gateaux 
differentiate at y. 

The existence of Gateaux derivative 5C[y; rj\ presupposes that: 

(i) C[y\ is defined; 

(ii) C[y -\- erj] is defined for all sufficiently small e. 
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Theorem 3.4. Let C,Gi, . . . , Qm be functionals defined in a neighborhood of y and 
having continuous Gateaux derivative in this neighborhood. Suppose that y is a weak 
local minimum of ([T]) subject to the boundary conditions ([2]) and the isoperimetric 
constrains Then, either: 
(ij^vj eC^a, j = l,...,m 

(4) 



SGi[y;vi] dgi[y;v2] 
Sg2[y;vi] Sg2[y;v2] 



5Q2 [y;f>i 



or 



SGmlmvi] 5Gm[y\v2] ■■■ 5Gm[y;vm 

(a) there exist constants Ai G M, « = 1, . . . , m for which 
(5) <5/:[y;77]-^A,<5g,[y;77] Vj; G C^, 



Proof. This proof is patterned after the proof of Troutman [101 Theorem 5.16]. Let 
us consider, for fixed directions 77, wi, W2, . . . , Wm, the auxiliary functions: 

l{p, qi,..., qrn) = C[y+pr] + qiVi H h 

gi(p, ^ Gi[y + PV + Qivi H \- qmVm], 



9miP, 91, ■ • • , Qm) = Gm[y + PV + 9l^'l H 1- qmVm], 

which are defined in some neighborhood of the origin in ]R"'+^, since C, Gi, 



■ , Gn 



themselves are defined in a neighborhood of y. Note that the partial derivative 

d d 
lp{p,qi, . . . = —l{p,qi, . . .,q,n) = —C[y + pri + qiVi H h q„iV„ ' 



dp 



lim^ 



C[y +{p + e)ri + qiVi -\ h g,„u,„] - C[y + prj + qiVi H h q,nV„ 



lim^-^i, 



C[y + er]] - C[y] 



with y = y+prj+qiVi- 
we have: 



Therefore, lp{p, qi,...,q 

IqAP^Hi ■ ■ ■ ^Ini) = SC[y;Vi], i 
{9j)p{P, 91: ■ ■ • , Qm) = SGj[y; v], 3 



i9j)q^ (p, 91: ■ ■ • , qm) = SGj [y; Vi\, i = 1, 



J 



: 1, 



d{l,gi, 



Hence, the Jacobian determinant 
is the following: 



evaluated at (p, qi , . 



5C[y; 77]. Similarly 
. . ,m, 

. . ,TO, 
. . . , TO. 

.,qm) = (0,0,. ..,0) 



(6) 



SC[y;i]] SCi[y;vi] 



5Ci[y;Vn 
5Gi [y;vr, 



SGmlmv] 5Gm[mvi] ■■■ SGm[y;vm 

Note also that the vector valued function {I, gi, . . . , gm) has continuous partial 
derivatives in a neighborhood of the origin, since C,Gi, ■ ■ ■ ,Gm have continuous 
Gateaux derivative in the neighborhood of y. With this preparation we can prove 
our theorem. Assume condition (i) does not hold for one set of directions: vi,V2, . . . ,Vr, 
and suppose there exists one direction rj for which the determinant © is nonva- 
nishing. Therefore, the classical inverse function theorem applies, i.e. the appli- 
cation {l,gi, . . . ,gm) maps a neighborhood of the origin in R™+^ onto a region 
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containing a full neighborhood of [C[y],Qi[y\, . . . ,Qm[y]). That is, one can find 
pre-image points (p, gi, . . . , and (p, gi, . . . , g„i) near the origin, for which the 
points y = y + pri + ^^iQiVi and y = y + pr/ + S™ j^^iUi satisfy the conditions: 

m<m<m, 

Gi[y] = Gt[y] ^ Gt[y], i = l,...,m. 

This shows that y cannot be a local extremal for C subject to constraints 
contradicting the hypothesis. Thus, for the specific set of directions; vi,V2, ■ ■ ■ ,Vrn 
the determinant ([6]) must vanish for each rj G C^^. We expand it by minors of the 
first column 

(7) dC[y;r]]-cof5C[y;r]]+Sgi[y;T]]-cof5gi[y;'q] + .. ■ + Sg,n[y;v]-cofSgm[y;v] =0, 

where we are using the notation cof to denote the cofactor. Dividing equation ^ 

Sgi[y;v.j] 



i, j = 1, . . . ,TO 



by cofSC[y; 77], since it is precisely the nonvanishing determinant 
and setting 

^ ^ cofSGiiy;!]] 
cofSC[y;'n] 

we obtain an equation equivalent to ((51). □ 



Note that condition (ii) of Theorem 13.41 can be written in the form 
(8) ^{^'f-t '^^) = ° ^ 

since the Gateaux derivative is a linear operation on the functionals (by the linear- 
ity of the ordinary derivative). 

Now, suppose that assumptions of Theorem 13.41 hold but condition (i) does not 
hold. Then, equation ([8]) is fulfilled for every 77 G C^^. Let us consider function 77 
such that r7(a) = 77(6) = and denote by J- the functional C — X]i=i ^iGi- Then we 
have 

= 5J'[y;rj] = —T[y + e7^]U^o 

= [\Fs{t,r{t),y''{t))Tj'^{t) + F,{t,r{t),y''{t))Tj'^{t))At, 



where the function F : [a,b]'' x M" x R" ^ R is defined by F{t, s, v) = L{t, s, v) 
J27Li KGi{t, s, v). Note that 



A 

F,(r,y-(r),y^(r))AT7y(i)) At = / F,(t, ^ (r), y^(T)Ar77(t)|*Z^ = 



/ a X'J a 

and 



s A 

F,iT,r{r),y^{T))ATrjit)] At 



a \J a 
b 



(^^ F,{r,r{r),y''{r))Ar^ rf (t) + F,(t, y-(r), ^^(t)) Ar77^ (<) | Ai 
rit),y^itWit) + f F,{r, r{r),y^{r))ATrj^{t)} At. 



Therefore, 

'F„(t,r(t),y^W)- / F,(T,r(r),y^(r))ATU^(t)At. 
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Since the function rj is arbitrary, Lemma 12.31 implies that 

J a 

for some c G and all t G [a, 5]^. Hence, 

(9) F^^{t,r{t),y''{t))^FS,r{t),y^{t)) 

for aU t e [a, b]''^ . 

We have just proved the following necessary optimality condition for the isoperi- 
metric problem with multiple constrains on time scales. 

Theorem 3.5. Let us assumptions of Theorem \3.4\ hold but condition ^ does not 
hold. If y E is a weak local minimum of the problem ([l])-([3]), then it satisfies 
the Euler- Lagrange equation Q for all t G [a, 6]*^ . 

3.2. Pareto optimality. Let us consider a finite number > 1 of (objective) 
functionals: 

(10) - / i.(i,y"(t),y^(i))Ai, z = l,...d, 

where a,b e T with a < b, L,{t,s,v) : [a,b]'' x R" x M" ^ R, i = 1, . . . d, have 
partial continuous derivatives with respect to the second and third variables for 
all t € [a,b]'', and Li{t,-,-) and theirs partial derivatives, i — l,...d, are rd- 
continuous at t. We would like to find a function y G C^^, satisfying the boundary 
conditions that renders the minimum value to each functional Ci, i = 1, . . . ,d, 
simultaneously. In general, there does not exist such a function, and one uses the 
concept of Pareto optimality. 

Definition 3.6. A function y G C^^ is called a weak locally Pareto optimal solution 
if there exists S > such that there does not exist y G C^^ with \\y — yWc^^ < ^ ^'i^d 

G {1, . . . ,d} : CM ^ CM A 3j G {1, . . . ,d} : CM < CM . 

Theorem 3.7. // y is a weak local minimum of the functional X^iLi li^i[y] with 
7i > for i ~ 1, . . . , d and X^iLi li ~ 1j then it is a weak locally Pareto optimal 
solution of the multiobjective problem with functionals (llOp . 

Proof. Let y be a weak local minimum of the functional X^iLi li^Av] with 7^ > 
for i = 1, . . . , d and X^iLi li ~ 1- Suppose on the contrary that y is not a weak 
locally Pareto optimal. Then, for every 6 > there exists y with \\y — yWc'-^ < ^ 
such that Vi G {l,...,d} we have CM] ^ CM] and 3j G {l,...,d} such that 
Cj[y] < Cj[y]. Since 7^ > for i = 1, . . . , d, we obtain J2i^i li^i{y] < SiLi li^M]- 
This contradicts our choice of y. □ 

Theorem 3.8. If y is a weak locally Pareto optimal solution of the multiobjective 
problem with functionals pO|) . then it minimizes each one of the scalar functionals 

CM, ie{l,...,d} 

subject to the constraints 

^j[y] ^ ^j[y], j = 1, • ■ ■ , and ji^i. 

Proof. Let ?} be a weak locally Pareto optimal solution of the problem on time scales 
(fTO|) and suppose the contrary, i.e. that for some i y does not solve the problem 
CM] ~* vnin subject to Cj\y] — Cj[y],j = 1, . . . ,d{j ^ i). Then, for every 5 > Q 
there exists y with \\y — y\\c^ < 5 such that Ci[y] < Ci[y] and Cj[y] = Cj[y],j = 
1, . . . ,d{j i). This contradicts the weak local Pareto optimality of y. □ 
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Example 3.9. Let T = {0,1,2}. We would like to find locally Pareto optimal 
solutions for 



Jo 

^2M= / {y{t + l)-2fAt 
Jo 



satisfying the boundary conditions y{0) = 0, y{2) — 0. Note that 

1 1 

t=0 t=0 

and that the possible solutions are of the form 



where a e M. On account of the above, we have Ci[y{t)] = a^, and £2 [2/(0] = 
4 + (a — 2)^. Using Theorem 13.71 we obtain that locally Pareto optimal solutions 
for functionals £1, £2 are 






if < 


= 


a 


if t 


= 1 





if t 


= 2, 







if t = 





yit) = 1 




if t = 


1, ae[0,2] 






if t = 


2 
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